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ABSTRACT 

We present the results of large-area 12 CO J— 3-2 emission mapping of three nearby field galax- 
ies, NGC 628, NGC 3521, and NGC 3627, completed at the James Clerk Maxwell Telescope as part 
of the Nearby Galaxies Legacy Survey. These galaxies all have moderate to strong 12 CO J=3-2 
detections over large areas of the fields observed by the survey, showing resolved structure and 
dynamics in their warm/dense molecular gas disks. All three galaxies were part of the Spitzer 
Infrared Nearby Galaxies Survey sample, and as such have excellent published multi-wavelength 
ancillary data. These data sets allow us to examine the star formation properties, gas content, 
and dynamics of these galaxies on sub-kiloparsec scales. We find that the global gas depletion 
times for dense/warm molecular gas in these galaxies is consistent with other results for nearby 
spiral galaxies, indicating this may be independent of galaxy properties such as structures, gas 
compositions, and environments. Similar to the results from the THINGS Hi survey, we do not 
see a correlation of the star formation efficiency with the gas surface density consistent with the 
Schmidt-Kennicutt law. Finally, we find that the star formation efficiency of the dense molecular 
gas traced by 12 CO J=3-2 is potentially flat or slightly declining as a function of molecular gas 
density, the 12 CO J=3-2/ J=l-0 ratio (in contrast to the correlation found in a previous study 
into the starburst galaxy M83), and the fraction of total gas in molecular form. 
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1. Introduction 

Gas and dust in disk galaxies are intimately 
linked to almost all aspects of their evolution, dy- 
namics, and appearance. For example, the in- 
terstellar medium (ISM) influences star formation 
and evolution, while these processes in turn have 
an effect on the distribution and turbulence of the 
ISM. One of the most important links is that be- 
tween the molecular clouds, which provide the fuel 
for star formation, and the location and rate of 
star formation. However the underlying physical 
process th at drives this relation is still p oorly un- 
derstood (jSchave fc Dalla Vecchial 120081). A re- 
cent detailed study by Lerov et al.l ( 2008 ). us- 
i ng the high resolu tion Hi maps from THINGS 
(j Walter et al.ll2008l ). found that many of the com- 
monly used prescriptions do not predict the star 
formation well. Tracing all components of the ISM 
across a wide variety of galaxies can lead to a 
greater understanding of the physics of the star 
formation process, its dependence on the prop- 
erties of the host galaxy's ISM, the evolution of 
galaxies, and the distribution of dark matter. 

One of the limiting factors in understanding the 
physics of star formation is our ability to trace 
the molecular gas clouds within galaxies. Much 
work has been done recently to examine the molec- 
ular gas within individual and small samples of 
nearby galaxies using various molecular transition 
lines that trace different gas densities and temper- 
atures. The most extensive observations thus far 



have used the millimeter 12 CO J— 1-0 transition 
to map molecular gas, including large surveys of 
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terferometers (e.g 


. Youne et al. 1995; Heifer et al. 


2003; Kuno et al. 


2007). This transition is most 



sensitive to cold, lower density molecular gas, 
which does not appe ar to trace star f ormation 



on a one-to-on e basis (IGreve et al.1 120051 ) . Recent 
studies such as Iono et al. ( 20091 ) have found that 



higher order transitions such as the 12 CO J=3- 
2 line more closely follow global star formation, 
nearly linearly over five orders of magnitude. 

The James Clerk Maxwell Telescope (JCMT) 
is a 15 meter single dish submillimeter telescope 
operated by the Joint Astronomy Centre (JAC) 
on Mauna Kea, Hawai'i, which after a recent ma- 
jor upgrade is undertaking the JCMT Legacy Sur- 
vejQ (JLS) a series of seven surveys from planetary 
to cosmological scales. One of these, the JCMT 
Nearby Galaxies Legacy Survey (NGLS), is the 
first large submillimeter survey of nearby galax- 
ies at ~15" spatial resolution. Taking advantage 
of new instrumentation on the JCMT, the survey 
will observe a well-selected sample of 155 nearby 
galaxies (within 25 Mpc) at 850 /im and 450 /Ltm 
and in the 12 CO J=3-2 line. The NGLS data set 
will be a powerful tool for studying the physics of 
the dusty ISM in galaxies as well as the interplay 
between star formation and the ISM. With a large, 
well-selected sample, we will be able to search for 
variations in the physical properties of the ISM as 
a function of galaxy type, metallicity, star forma- 
tion rate, mass, and environment. By having lim- 
ited our sample galaxies to distances closer than 
25 Mpc, we will be able to study spatial scales 
of 0.2-2 kpc and to search for variations inside a 
single galaxy as well as between galaxies. 

Since November 2007 the first phase of the sur- 
vey has been complete, observing 57 galaxies of 
the sample in the 12 CO J=3-2 line; 21 overlap 
with the Spi tzer Infrared Nearby Galaxies Sur- 
vey fSINGS. iKennicutt et aill2003l) sample, and 
36 are members of the Virgo Cluster. 12 CO J=3- 
2 observations for the remaining SINGS and field 
sample g alaxies are cu r rently underway. Following 
on from IWilson et aD ( 20091 hereafter Paper I), 



which investigated four SINGS subsample galaxies 
that are members of the Virgo Cluster, we present 



t http: / / www. j ach. hawaii . edu /JCMT/ surveys / 
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Table 1 

Summary of Previously Measured Galaxy Properties. 



Name 


a(J2000.0) 


5(J2000.0) 


Type 




d 


M B 


Mmi 


Mm 2 


log L Ha 


i 










(kms -1 ) 


(Mpc) 


(mag) 


(10 8 M B ) 


(10 s M@) 


(log org s -1 ) 


n 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


en) 


NGC628 


01 36 41.77 


+ 15 4700.5 


SA(s)c 


648 


7.3 


-19.58 


38.0 


6.3 


40.87 


7.0 


NGC 3521 


1105 48.88 


-00 02 05.7 


SAB(rs)bc 


778 


10.7 


-19.85 


80.2 


19.5 


41.00 


72.7 


NGC3627 


1120 15.03 


+ 12 59 29.6 


SAB(s)b 


697 


9.4 


-20.44 


8.18 


41.2 


41.11 


61.8 



Note. — Col. (1): Galaxy name. Col. (2) and (3): J2000.0 right ascension and declination as given in RC3 (dc Vaucoulcurs ct al. 1991). Units 
of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutcs, and ar cseconds. Col. (4): Galaxy morphology 
type from NED. Col. (5): S ystemic velocity from the Hi line (Heliocentric). Col. (6): Ad opted distance { jKarachcntscv ct al. 2004; Walte r" et al.1 
2008; Frccdman ct al. 2001, respectively). Col. (7): Total absolute B-band magnitude JKennicutt et all 120081) . Col. (8): The H I mass from 
THIN GS Walter et alj|2008fl . Col. (9): T he H 2 mass from the N RAO 12 m 12 CO J=l-0 on-the-fly obs ervations for BIMA SONG fHelfer et all 
I2003D . Col. (10): The Ha luminosity from lKennicutt et aD ||2008I 1. Col. (11): The adopted inclination llWalter et al.ll2008l : Ide Blok et al.H2008T) 



the 12 CO J=3-2 results for three SINGS galaxies 
outside the cluster environment. 

For our initial examination of the 12 CO J=3- 
2 data in non-cluster NGLS galaxies, we have 
chosen to look at three well known local spi- 
ral galaxies that were completed in the early 
stages of the survey: NGC 628, NGC 3521, and 
NGC 3627. All three of these galaxies have ex- 
tensive multi-wavelength ancillary data, most re- 
cently from SINGS and the follow-up surveys as- 
sociated with it. Having been included in SINGS, 
all three have publicly available infrared obser- 
vations (seven bands from 3.6 to 160 /an, plus 
some spectroscopy). All three were also included 
in The Hi Nearb y Galaxy Survey (THINGS, 
Walter et ai1l2008l) . the VLA Hi follow up sur- 
vey to SINGS. Additionally, they were all ob- 
served in the 12 CO J=l-0 line on the Berkeley- 
Illinois-Maryland Association (BIMA) millime- 
ter interferometer as part of the B IMA Survey 
of N earby Galaxies (BIMA SONG, iHelfer et al 



2003h . and two (NGC 3521 and NGC 3627) have 
single dish 12 CO J— 1-0 observations taken with 
the N obeyama 45-m telescope by iKuno et~al 
(120071) . Comparing some of the abundant an- 



cillary data available with our JCMT 12 CO J=3-2 
observations provide us with a powerful tool for 
exploring star formation within these galaxies. 
Table Q] outlines some of the previously published 
properties of these three galaxies. 

NGC 628 (M 74) is an isolated, near face-on spi- 
ral of type SA(s)c. It is known to have an extended 
Hi disk, extending to appro ximately three times 
the optical Holmberg radius ([Kamphuis fc Briggs 



19921 ). 12 CO J =l-0 observations of NGC 628 for 



BIMA SONG (IHelfer et all 120031 ) only cover the 
inner region of the galaxy, but show that the 
molecular gas is clumpy and distributed mostly 
along the spiral arms. It is slightly less luminous 
than the other two galaxies we are examining here, 
and also has the lowest Ha and 12 CO J=l-0 lumi- 
nosity of the three. We ado pt the same distanc e 
estimates in this paper as in lWalter et al.l ([20081 ). 
For NG C 628 the distance w e use i s 7.3 Mpc, de- 
rived by iKarachentsev et al.l (|2004l ) using the lu- 
minosity of the brightest stars. 

NGC 3521 is a highly inclined flocculent spiral 
galaxy (type SAB(rs)b c) that shows a r ing-like 
i 2 CO J= 1-0 structure ([Heifer et al.ll2003l) . It is 
associated with the nearby Leo Triplet (of which 
another of our galaxies, NGC 3627, is a member). 
It has a moderatel y high Hi mass (8 x 10 9 A4q, 
Walter et all 120081) . the highest of these three 
galaxies, and the H I disk extends well beyond the 
stellar disk. The adopted dista nce to NGC 3521 
is 10.7 Mpc (|Walter et al.ll2008l Hubble flow dis- 
tance using the velocity listed on NED). 

NGC 3627 (M 66) is an active, asymmetric 
barred spiral galaxy (type SAB(s)b), a mem- 
ber of the Leo Triplet well known for its un- 
usual kinematics that are infl uenced by both its 
bar and external i nteractions jzhang et al] 19931: 
Reuter et al.lll996l : Regan et al.l l2002: Ch emin et al 



2003). It is the most luminous and most ac- 
tive of these three galaxies, and has the high- 
est molecular gas content of the three as de- 
duced from 12 COJ= 1-0 observations (A^Ho = 
(4.1 ±0.4) x 10 9 M e , IHelfer et al.ll2003h . In con- 
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trast, the H I mass is the lowest of the three galax- 
ies, such that the molecular gas dominates the gas 
content of the galaxy. The ado pted distance to 
NGC 3627 is 9.3 Mpc, derived bv lFreedman et al 
(|2001l ) using Cepheid variables. 



In this paper we focus on the molecular gas 
kinematics and detailed star formation properties 
of these three galaxies. § [2] describes our obser- 
vations at the JCMT and data reduction process, 
and the processing of ancillary data. § [3] presents 
the resulting moment maps from our 12 CO J=3-2 
observations. § |4] compares our new 12 CO J=3-2 
observations to existing Hi and 12 CO J=l-0 data. 
§[5] deals with the dynamics of the molecular gas, 
looking at differences in the velocity fields from 
molecular and atomic gas, and comparing the ro- 
tation curves we derive from our 12 CO J=3-2 ob- 
servations with thos e from the H I data of THINGS 
( Walter et al. 20081 ). In §[6j we discuss several as- 
pects of star formation within these galaxies. Fi- 
nally we present our conclusions in § [71 



2. Observations and Data Reduction 
Overview 

2.1. JCMT 12 COJ=3-2 Data 

The HARP-B 12 CO J=3-2 line (rest frequency 
345.79 GHz) observations for the three galaxies 
we are examining here, NGC628, NGC3521, and 
NGC 3627, were taken over multiple runs between 
November 2007 and March 2008 as part of the 
JLS program. For all three galaxies we mapped 
a rectangular area corresponding to -D25/2 us- 
ing a basket weave raster scanning method (see 
Appendix [A] for details). We set up the back- 
end spectrometer, the Auto-Correlation Spec- 
trometer Imaging System (ACSIS), with a band- 
width of ~1 GHz and a resolution of 488.28 kHz 
(0.423 kms -1 at the frequency of the observa- 
tions). In each scan, we integrated for 10 seconds 
per pointing within the target field (if all 16 recep- 
tors were in working order, see Appendix [A]) , and 
scans were repeated until we reached the target 
RMS for the survey 19 mK (T^) after binning 
to 20 kms -1 resolution. Individual setup and ob- 
serving details for the three galaxies are listed in 
Tabled 

Data reduction and most of the analysis were 
carried out with recent versions (Humu and 



Lehuakona) of the Starling software package that 
is maintained by the J AC. Additionally, analy- 
sis tasks in MIRIAD were used after the creation 
of cubes and moment maps. The general reduc- 
tion procedure was to inspect the raw data files, 
flag poor data where necessary, combine the raw 
data into a cube, subtract a baseline from the 
cube, and collapse it into moment maps. We con- 
verted the 12 CO J=3-2 data from T^ to the main 
beam temperature scale by dividing by T)mb = 0.6. 
Full details of the observations and data reduction 
methods are described in Appendix El 

2.2. Ancillary Data 

All three galaxies in this paper are well studied 
members of the SINGS sample, and as a result 
have a large quantity of published ancillary data, 
from the far ultraviolet to centimeter radio ob- 
servations. We obtained 24 /x m and Ha maps 
from the SINGS survey itself dKennicutt et al 



2003) at the SINGS websited. Hi data cubes and 
0th moment maps came from the THINGS sur- 
vey ((Walter et all l2008t ) websiteH There were 
two sources for 12 CO J=l-0 maps , single dish 
Nobeyama maps from iKuno et al. |2007jl for 
NGC 3521 and NGC 3627, and the BIMA SONgEI 
12 CO J=l-0 maps for NGC 628 as there were no 
Nobeyama data available. Although BIMA SONG 
data was available for all three, they only have 
good sensitivity in the central part of the field 
covered by the NGLS, so we preferred the single 
dish Nobeyama data where available as they cover 
a similar field size, and the Nobeyama beam size 
for 12 CO J=l-0 data is well matched to the JCMT 
at 12 CO J=3-2 (so no convolution is required). In 
all cases, the available ancillary data maps were 
of similar or higher resolution than the 14'.'5 beam 
of the NGLS 12 CO J=3-2 observations, and they 
cover a field-of-view of similar size or larger. 

The procedures for processing the SINGS 
24 /im and Ha data are outlined in Paper I, and 
follows the notes and conversions supplied the 



available for download from 

http: / /starlink.jach. hawaii.edu/ 
§ http: / /sings. stsci.edu 

^http: / /www. mpia-hd.mpg.de/THINGS/Overview. html 
" http: / / www.nro.nao.ac.jp/ nro45mrt/COatlas / 
"Obtained from NED, http://nedwww.ipac.caltech.edu/ 
level5/March02/SONG/SONG.html 
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Table 2 

Summary of Observing Parameters. 



Name 

(i) 


Obs. Dates 
(2) 


PA 

(deg) 
(3) 


D 2 5, ma j/2 
(arcscc) 
(4) 


D 25 , min /2 
(arcscc) 
(5) 


T sys 
(K) 
(6) 


r 

(225 GHz) 
(7) 


w 

(8) 


AT 

(mK) 
(9) 


Ref. Pos. 
(arcscc) 
(10) 


NGC628 


071123, 071124, 071125 
080107, 080113 


25 


312 


288 


388 


0.097 


80 


15.8 


432 


NGC3521 


080310, 080311 


163 


324 


162 


402 


0.109 


65 


16.8 


444 


NGC3627 


080112, 080113 


173 


264 


120 


328 


0.084 


70 


13.2 


384 



Note.— Col. (1): Galaxy name. Col. (2): Date of observations in 2007 and 2008, YYMMDD. Col. (3): Position angle of the major 
axis from the optical component, and orientation of the observed field. Col. (4): Major axis length of the target scan region, D25,maj /2. 
Col. (5): Minor axis length of the target scan region, D25,min f2- Col. (6): Median system temperature over all observations. Col. 
(7): Median atmospheric opacity at 225 GHz as measured with the JCMT water vapor radiometer. Col. (8): Typical integration time 
per point in the map. Col. (9): Pons noise (Tj£) in line-free regions of the spectra at 20 krns" 1 resolution. Col. (10): Location of the 
shared off position, in arcseconds cast of the center. 



SINGS User's Guide lSINOS team! (|2007t) . To com- 
pare these data to the NGLS 12 CO J=3-2 maps, 
we first had to convolve them to give a smoothed 
image with a 14'.'5 circular beam to match the 
JCMT resolution at 345 GHz. The Ha images 
contained a significant (and uneven) background 
sky level that needed to be subtracted by fitting a 
surface to the image. The fit was limited to a first 
order polynomial as any higher orders would be- 
gin to trace low level Ha emission from the galaxy. 
Several artifacts caused by over-exposed stars were 
also removed from these images using upper and 
lower thresholds prior to smoothing. The 24 fim 
data did not require any special processing prior 
to convolving the map to the JCMT resolution. 
The same smoothing procedure was performed on 
the BIMA SONG 12 CO J=l-0 map for NGC628; 
however, the Nobeyama 12 CO J— 1-0 data for the 
other two galaxies were not convolved to this beam 
as the angular resolution already closely matches 
the JCMT. 

For the H I data, we use the naturally weighted 
zeroth moment maps that are available on the 
THINGS Public Data Repository web site 



(http://www.mpia-hd.mpg.de/THINGS/Data.html) 
for comparison to our 12 CO J=3-2 data. These H I 
data are well matched in resolution to the NGLS 
observations, with our single dish 12 CO J=3-2 ob- 
servations being at only slightly poorer resolution 
than the naturally weighted interferometric H I 
data (14'/5 for the JCMT compared to 9" to 14" 
for the VLA data). As with the other ancillary 



data, we convolved the H I intensity (Oth moment) 
maps to match the JCMT resolution. 

Unlike the Hi intensity maps, we cannot di- 
rectly convolve the beams of the velocity fields 
(1st moment) and velocity dispersion (2nd mo- 
ment) maps from the THINGS survey for compari- 
son to our 12 CO J=3-2 moment maps. So instead, 
we obtained the naturally weighted THINGS Hi 
data cubes for the three galaxies, convolved the 
cubes to a 14'/5 circular beam, and then created 
moment maps from the resulting cubes (using a 
similar procedure to the creation of the 12 CO J=3- 
2 maps). While this also produced Oth moment 
maps, we prefer to use the convolved version of 
the THINGS maps in this case as their moment 
map production procedure has been refined for 
their data set. After convolving to the appropri- 
ate resolution, all ancillary data were aligned to 
the JCMT 12 CO J=3-2 map and resampled to the 
same pixel scale as the 12 CO J=3-2 maps (~7'/3 
pixels). 

3. Results: 12 CO J=3-2 Images 

The 12 COJ=3-2 moment maps we obtained 
from our HARP-B raster scan mapping are shown 
in Figures ED (NGC628), EH (NGC3521), and 
(NGC3627). Each figure contains four panels, 
showing the integrated 12 CO J=3-2 intensity over- 
laid on the optical DSS II (R band) image, the inte- 
grated 12 CO J=3-2 intensity by itself, the velocity 
field, and the velocity dispersion map. The sec- 
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Fig. 1. — 12 CO J=3-2 moment maps of NGC628. a) Integrated 12 CO J=3-2 intensity distribution (moment 
0) contours for NGC628 overlaid onto an optical DSSn R band image. Contour levels are 0.5, 1, and 
2 Kkms -1 (temperatures in all moment maps are T m b). b) Integrated 12 CO J=3-2 intensity contours as 
in first panel but with our 12 CO J=3-2 map in the background instead of the DSS II image. White is low (or 
blank), and black is high. The blue box in this panel shows our target mapping region (note that the useful 
region of the final image is slightly larger than these boxes, see Appendix [A"|l . c) 12 CO J=3-2 velocity field 
(moment 1) for NGC628. Contour levels (from blue to red end) are 636, 646, 656 (thick contour, systemic 
velocity), 666, and 676 kms -1 . d) 12 CO J=3-2 velocity dispersion map (moment 2) for NGC628. Contour 
levels are 4, and 6 kms -1 . 
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Fig. 2. — 12 CO J— 3-2 moment maps of NGC 3521. a) Integrated 12 CO J=3-2 intensity distribution (moment 
0) contours for NGC 3521 overlaid onto an optical DSSn R band image. Contour levels are 0.5, 1, 2, 4, and 
8 Kkms -1 . b) Integrated 12 CO J=3-2 intensity contours as in first panel but with our 12 CO J=3-2 map in 
the background instead of the DSS n image. White is low (or blank), and black is high. The blue box in this 
panel shows our target mapping region, c) 12 CO J=3-2 velocity field (moment 1) for NGC 3521. Contour 
levels (from blue to red end) are 608, 658, 708, 758, 808 (thick contour, systemic velocity), 858, 908, 958, 
and 1008 kms -1 . d) 12 CO J=3-2 velocity dispersion map (moment 2) for NGC 3521. Contour levels are 5, 
10, 20, and 30 kms" 1 . 
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Fig. 3. — 12 C0 J— 3-2 moment maps of NGC 3627. a) Integrated 12 C0 J=3-2 intensity distribution (moment 
0) contours for NGC 3627 overlaid onto an optical DSS n R band image. Contour levels are 0.5, 1, 2, 4, 8, 16, 
32, and 64 Kkms -1 . b) Integrated 12 C0 J=3-2 intensity contours as in first panel but with our 12 CO J=3-2 
map in the background instead of the DSSn image. White is low (or blank), and black is high. The blue 
box in this panel shows our target mapping region, c) 12 CO J=3-2 velocity field (moment 1) for NGC 3627. 
Contour levels (from blue to red end) are 555, 595, 635, 675, 715 (thick contour, systemic velocity), 755, 795, 
835, and 875 kms -1 . d) 12 CO J=3-2 velocity dispersion map (moment 2) for NGC 3627. Contour levels are 
4, 8, 16, 32, and 64 kms" 1 . 
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ond panel of each figure shows the target field for 
the survey in blue, a D 2 5,maj/2 by D 2 5,mm/2 box, 
where theoretically we have complete sensitivity 
coverage to our target depth. In practice, signal- 
to-noise variations between receptors, and miss- 
ing/flagged receptors in the array, make the sen- 
sitivity coverage of the field more complex. Addi- 
tionally, the region with good sensitivity is slightly 
larger than the target field due to the over-scanned 
area (see Appendix |A"|) . At the adopted distances 
of the three galaxies, the 14'.'5 beam of the JCMT 
12 CO J=3-2 observations gives us resolutions of 
510, 750, and 660 pc for NGC 628, NGC 3521, and 
NGC 3627, respectively. 

The integrated intensity maps of NGC 628 
(Fig. QJ,,b) show that the warm molecular gas 
detected in the NGLS 12 CO J=3-2 line observa- 
tions is patchy, of relatively low intensity, and 
mostly follows the spiral arms of the galaxy. The 
center region shows neither a strong concentra- 
tion nor a strong dip in the molecular gas. The 
few higher peaks in the map are mostly dis- 
tributed along the spiral arms, the strongest at 
3.7 ± 0.7 Kkms -1 being ~1.'5 east of the galaxy 
center (though there are several others of similar 
strength). The total 12 CO J=3-2 luminosity for 
NGC 628 is (3.1 ± 1.3) x 10 7 Kkms" 1 pc 2 . 

NGC 6 28 is almost face on (inclination of 
about 7°. IWalter et~aT1l2008h . so the velocity field 
(Fig.nj) shows only a narrow velocity range across 
the galaxy (^50 kms -1 ). As our target field only 
covers the center of the galaxy, the coverage is 
patchy, and the inclination is so low, we can- 
not tell from the 12 CO J=3-2 data alone if we 
have data far enough out to reach the flat por- 
tion of the rotation curve, though as we only map 
out to D25/2 that is unlikely. Fig. [TJi shows the 
12 CO J=3-2 velocity dispersion map (moment 2 
map) for NGC 628. The velocity dispersion is 
generally low, around 2 to 6 kms -1 , with the 
higher velocity dispersion around the integrated 
12 CO J=3-2 intensity peaks. The median veloc- 
ity dispersion is 3.2 kms -1 . These velocity dis- 
persions are similar to the internal velocity dis- 
persions seen in gia nt molecular clouds (GMC) 
in the Milky Way (|Solomon et all Il987t ). This 
suggests that the main contribution to velocity 
dispersion in the molecular gas for NGC 628 is 
internal dispersion within GMCs, rather than 
inter-cloud velocity differences. Alternatively, 



ICombes fc BecauaertJ (1997) also looked at the 
velocity dispersion of molecular gas in this galaxy, 
using IRAM 30 m telescope observations at 
12 CO J=l-0 and 12 CO J=2-l (23" and 12" res- 
olution, respectively), and finding a near con- 
stant velocity dispersion as a function of radius 
of ^6 kms -1 they conclude that cloud-cloud ve- 
locity dispersions dominate over internal disper- 
sions. Without resolving the clouds themselves 
however we cannot tell for certain which compo- 
nents contribute most to the velocity dispersion 
measured within our beam. We would expect 
contributions from the rotation velocity gradient 
across our beam to be minimal beyond the center 
of the galaxy given the low inclination. 

The integrated 12 CO J=3-2 intensity maps for 
NGC 3521 (Fig.[2h,b) show that its molecular gas 
is concentrated towards the inner regions of the 
optical disk, with a slight central dip forming a 
torus- like structure in the inner ~1' (~3 kpc along 
the major axis). The high inclination angle of the 
galaxy makes it difficult to identify other struc- 
tures within the molecular gas disk. The detected 
12 CO J=3-2 disk extends along the major axis out 
to a radius of ~130", or ^6.7 kpc, which is about 
one third of the optical radius I?25,min- There are 
two peaks in the map on either side of the central 
dip (<30" from the center along the major axis 
in both directions) that are of about equal inten- 
sity in our 12 CO J=3-2 data, with the northern 
peak at 14.6 ± 0.6 Kkms -1 and the southern at 
15.1 ±0.7 Kkms -1 . The total 12 CO J=3-2 lumi- 
nosity for NGC 3521 is (2.8 ± 0.3) x 10 s Kkms -1 
pc . 

The 12 COJ=3-2 velocity field for NGC 3521 
(Fig. [5]:) shows the general characteristics of a 
moderate to highly inclined rotating disk galaxy 
(the well known spider diagram pattern). Over 
the region where we detect 12 CO J=3-2 there do 
not appear to be any significant distortions in the 
velocity field. A brief examination of position- 
velocity slices through the cube for NGC 3521, 
along and parallel to the major axis, showed 
no significant detections of molecular gas with 
anomalous velocities, and that the rotation curve 
appears to flatten out within the radius of the 
disk detected in 12 CO J=3-2. The former sug- 
gests that the 12 CO J=3-2 velocity field is tracing 
the bulk motions of the gas disk. Fig. [2ji shows 
the 12 CO J=3-2 velocity dispersion map (moment 
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2 map) for NGC3521. It has a distinctive 'X' 
pattern through the center of the galaxy (best 
seen with the 20 kins -1 contour), which seems to 
be aligned with the velocity field contours. The 
highest velocity dispersions are seen towards the 
center of the galaxy, where the velocity field con- 
tours are spaced closest together. Given the high 
inclination of this galaxy, these may indicate that 
the velocity dispersion map is partially indicating 
regions where there is a large component of the 
rotational velocity gradient in the direction of the 
line of sight. The median velocity dispersion is 
15.0 kms" 1 . 

The integrated intensity maps of NGC 3627 
(Fig. [3^, b) shows a distinct pattern of an inclined 
barred spiral. Molecular gas is found throughout 
much of the optical disk, with 12 CO J=3-2 de- 
tected right up to and sightly beyond the edges 
of our target field (the blue box in the upper right 
panel). It is mostly constrained to the bar and 
along the spiral arms, with no gas detected in the 
intervening regions. The spiral arms are asymmet- 
rical in 12 CO J=3-2, with the arm that extends to 
the south from the western side of the galaxy at 
least 50% longer than the arm that extends north 
from the eastern side. The center and the ends 
of the bar have particularly strong concentrations 
of 12 CO J=3-2, with the peak in the bar/galaxy 
center at 81 ± 1 Kkms -1 . The peak intensity 
at the southern bar end is 43.0 ± 0.6 Kkms -1 , 
while in the northern bar end the peak is at 
30.2 ± 0.6 Kkms" 1 . The total 12 CO J=3-2 lumi- 
nosity for NGC 3627 is (3.1 ±0.2) x 10 s Kkms" 1 
pc 2 . 

Fig. [3t shows the 12 CO J=3-2 velocity field 
for NGC 3627. The field shows some distortion 
between the inner bar and spiral arms, with the 
gradient of the field in the inner bar region aligned 
close to north-south (the receding side is ^180° 
East of North), while further out the alignment 
is closer to ~170°. This resembles s treaming mo- 



tions associated wi th spiral arms ([Bosmal 11981 
Knapen et al. 19931) . Additionally, the kinemat 



ical major and minor axes are not orthogonal, 
characteristi c of oval dist ortions (such as bars) 
or warping ( Bosmal Il981 ). Fig. [2ji shows the 
12 CO J=3-2 velocity dispersion map (moment 
2 map) for NGC 3627. The velocity dispersion 
is moderate to low in the spiral arms (~7 to 
20 kms -1 ), but much higher in the bar region, 



particularly in the center of the galaxy where the 
velocity dispersion is over 80 kms -1 , though there 
is a large velocity gradient across the 14'.'5 JCMT 
beam in this region. The median velocity disper- 
sion is 13.0 kins -1 . 



4. 



The 12 COJ=3-2/J=l-0 Ratio and the 
Distributions of Molecular and Atomic 
Gas 



The 12 CO J=3-2/J=l-0 ratio is an important 
indicator of the portion of warm/dense molecu- 
lar gas relative to the total molecular gas. As 
noted in Paper I it is also important in the calcu- 
lation of the molecular gas mass, as to date con- 
versions from CO luminosity have been done with 
the 12 CO J=l-0 line. We calculated 12 CO J=3- 
2/J=l-0 maps for all three galaxies by the same 
method as Paper I, taking the NGLS 12 CO J=3-2 
integrated intensity maps (in main beam temper- 
ature) and dividing by the available 12 CO J=l- 
maps of the same resolution (Nobeyama maps 
for NGC 3521 and NGC 3627, and the smoothed 
BIMA SONG 12 COJ=1-0 maps for NGC 628). 
The resulting 12 CO J=3-2/ J=l-0 ratio maps are 
shown in the upper left panels of Figs.|U[5j and El 
for NGC 628, NGC 3521, and NGC 3627, respec- 
tively. 

The 12 CO J=3-2 detection in NGC 628 is rel- 
atively weak compared to the other two sources, 
and the rapid drop off in the BIMA SONG map- 
ping sensitivity beyond ~lf5 (^3.2 kpc) from 
the galaxy center reduces the useful area of the 
12 CO J=3-2/J=l-0 ratio map. Within the use- 
ful region, the gas distribution traced by both 
12 CO transitions follow each other closely. The 
12 CO J— 3-2/ J =1-0 ratio is generally low (mostly 
between 0.15 and 0.45) with the exception of a few 
peaks (some of which may be due to poor sensi- 
tivity in one tracer or the other). It is difficult 
to discern a trend with the limited coverage we 
have for this galaxy. For NGC 3521 the 12 CO J=3- 
2/ J =1-0 ratio map shows that throughout most of 
the galaxy the 12 CO J=3-2/J=l-0 is in the range 
0.25 to 0.5. In the center of the galaxy, where the 
NGLS 12 CO J=3 -2 map shows a st rong dip (see 
Fig. E}, and the iKuno et all (|2007h 12 CO J=l-0 
also shows a central depression, the ratio drops be- 
low 0.25. Just to the southwest of the center the 
ratio peaks over 0.75, with some stronger peaks 
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Fig. 4. — Comparisons of the NGLS 12 CO J— 3-2 data to other atomic and molecular gas tracers for NGC 628. 

a) The ratio of the NGLS integrated 12 CO J=3-2 intensity to the BIMA SONG integrated 12 CO J=l-0 in- 
tensity. The magenta, black, and cyan contours mark where the ratio is 0.15, 0.30, and 0.45, respectively. 

b) The H2 (CO J=3-2) molecular gas mass from the NGLS observations as a fraction of the total gas mass 
(H2 [CO J=l-0] mass plus Hi mass from THINGS). The magenta, black, and cyan contours mark where H2 
(CO J=3-2) makes up 25%, 50%, and 75% of the total gas mass, respectively, c) The velocity field differ- 
ences between the THINGS Hi maps and the NGLS 12 CO J=3-2 maps (Hi velocity field minus 12 CO J=3-2 
velocity field). Redder values represent positive differences (Hi velocity greater than the 12 CO J=3-2 veloc- 
ity), while bluer values represent negative differences (Hi velocity less than the 12 CO J=3-2 velocity). The 
cyan, black, and orange contours represent differences of +6, and —6 kms -1 , respectively, d) The ratio 
of the NGLS 12 CO J=3-2 velocity dispersion to the THINGS Hi velocity dispersion. The magenta, black, 
and cyan contours mark where the 12 CO J=3-2 velocity dispersion is 25%, 50%, and 75% of the Hi velocity 
dispersion, respectively. 
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Fig. 5. — Comparisons of the NGLS 12 CO J=3-2 data to other atomic and molecular ga s trace rs for 



(2007) inte- 



NGC3521. a) The ratio of the NGLS integrated 12 C0 J=3-2 intensity to the Kuno et al. 
grated 12 CO J=l-0 intensity. The magenta, black, and cyan contours mark where the ratio is 0.25, 0.50, 
and 0.75, respectively, b) The H2 (CO J=3-2) molecular gas mass from the NGLS observations as a fraction 
of the total gas mass (H2 [CO J=l-0] mass plus Hi mass from THINGS). The magenta, black, and cyan 
contours mark where H 2 (CO J=3-2) makes up 25%, 50%, and 75% of the total gas mass, respectively, c) 
The velocity field differences between the THINGS Hi maps and the NGLS 12 CO J=3-2 maps (Hi velocity 
field minus 12 CO J=3-2 velocity field). Redder values represent positive differences (Hi velocity greater 
than the 12 CO J=3-2 velocity), while bluer values represent negative differences (Hi velocity less than the 
12 CO J=3-2 velocity). The cyan, black, and orange contours represent differences of +40, and —40 kms -1 , 
respectively, d) The ratio of the NGLS 12 CO J=3-2 velocity dispersion to the THINGS Hi velocity disper- 
sion. The magenta, black, and cyan contours mark -vfbere the 12 CO J=3-2 velocity dispersion is 20%, 40%, 
and 60% of the H 1 velocity dispersion, respectively. 
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Fig. 6. — Comparisons of the NGLS 12 CO J=3-2 data to other atomic and molecular ga s trace rs for 



NGC3627. a) The ratio of the NGLS integrated 12 CO J=3-2 intensity to the Kuno et al.l (|2007t l inte- 
grated 12 CO J =1-0 intensity. The magenta, black, and cyan contours mark where the ratio is 0.15, 0.30, 
and 0.45, respectively, b) The H2 (CO J=3-2) molecular gas mass from the NGLS observations as a fraction 
of the total gas mass (H2 [CO J=l-0] mass plus Hi mass from THINGS). The magenta, black, and cyan 
contours mark where H 2 (CO J=3-2) makes up 25%, 50%, and 75% of the total gas mass, respectively, c) 
The velocity field differences between the THINGS Hi maps and the NGLS 12 CO J=3-2 maps (Hi velocity 
field minus 12 CO J=3-2 velocity field). Redder values represent positive differences (Hi velocity greater 
than the 12 CO J=3-2 velocity), while bluer values represent negative differences (Hi velocity less than the 
12 CO J=3-2 velocity). The cyan, black, and orange contours represent differences of +30, and —30 kms -1 , 
respectively, d) The ratio of the NGLS 12 CO J=3-2 velocity dispersion to the THINGS Hi velocity disper- 
sion. The magenta, black, cyan, green, and pale blue fgntours mark where the H2 velocity dispersion is 25%, 
50%, 75%, 100%, and 150% of the Hi velocity dispersion, respectively. 



around the edges of the map. For NGC 3627, the 
general 12 CO 7=3-2 structure matches we l l with 



both the BIMA SONG and iKuno et all (|2007h 
12 CO J— 1-0 maps of the galaxy, though there is 
more inter-arm gas in those maps than in the 
NGLS map. The 12 CO J=3-2/ J=l-0 ratio is low 
throughout most of the galaxy where we detect 
12 CO 7=3-2, only going above 0.45 in four loca- 
tions (the center, the bar ends, and the concentra- 
tion to the south of the southern bar end). The 
median 12 CO J=3-2/ J =1-0 ratios for NGC628, 
NGC 3521, and NGC 3627 are 0.35, 0.33, and 
0.18, respectively (for pixels where the 12 CO 7=3- 
2 signal-to- noise is 3 or greater). 

The standard method to calculate the column 
density of H2 gas using 12 CO 7=1-0 as a tracer is 
to adopt a CO-to-H 2 conversion factor (Xco)- As 
in Paper I we ar e adopting a value o f 2 x 10 2 



(K km s^ 1 )" 1 (|Strong et alJ ll988). In order to 



calculate the molecular gas mass from 12 CO J=3 
2 measurements we need to adopt a 12 CO J=3 
2/J=l-0 value to change our 12 CO J=l-0-to-H 2 
conversion factor in to a 12 CO 7=3-2-to-H 2 con- 
version factor. The choice of this conversion ratio 
depends strongly on where we assume the molecu- 
lar gas that we are tracing with 12 CO J=3-2 is lo- 
cated. Our assumption is that 12 CO 7=3-2 traces 
pre-star formation molecular gas that is in gen- 
eral denser than 12 CO J=l-0. In Paper I we used 
a constant ratio of 0.6, representative of ratios 
seen in Galactic and extragalactic GMCs (based 
on the results of iLerov et alj|2005 Wilson et al 



Il999t llsraell l2009al )bl). essentially focusing on the 
molecular gas closest to the star formation (un- 
like the above median 12 CO J=3-2/ J=l-0 ratios 
for the three galaxies, which are more representa- 
tive of the global molecular gas population). The 
limitation of this assumption is that the cause 
of the differences between different molecular gas 
tracers is degenerate, thus higher 12 CO 7=3-2 val- 
ues could indi cate warmer rat her than denser gas 
(as found by lOka et al.l 20071) . It has also been 



suggested (jTosaki et alj 120071) some 12 CO 7=3-2 



emission may be due to post-star formation heat- 
ing of the molecular gas. In this paper, where 
we calculate an H2 mass using this assumption 
( 12 CO J=3-2/J=l-0 = 0.6) we will refer to 'H 2 
(CO J=3-2).' 

If we instead wish to examine the more dif- 
fuse general molecular gas content in a galaxy 



then the available 12 CO 7=1-0 data is a more ap- 
propriate tracer. So in cases where we require 
the total gas content we use this instead, using 
the adopted Xco value to obtain the H 2 gas col- 
umn density, and refer to the 'H 2 (CO J=l-0)' gas 
mass. As the only 12 CO J=l-0 observations we 
have for NGC 628 come from BIMA SONG, and 
given the limitations of those data described above 
(see §2.2p . we have only used parameters needing 
H 2 (CO 7=1-0) for the other two galaxies in our 
later analysis. 

In order to compare the dense molecular gas 
content to the total gas (Hi plus H 2 [CO 7=1- 
0]), we produced maps of the H 2 (CO J=3-2) gas 
mass surface density fraction (Sh 2 (co j=3-2) di- 
vided by £hi+h 2 (co j=i-o))j where the H 2 surface 
densities were derived as above, and the H 1 surface 
density came from the THINGS integrated H 1 in- 
tensity maps smoothed to the JCMT beam. The 
resulting H 2 gas mass fraction maps are shown in 
the upper right panel of Figs. [4] [5j and [6l Since 
^h 2 (coj=3-2) is n °t in the denominator for this 
calculation it should be noted here that there are 
some pixels in these maps where this fraction is 
greater than 1. 

As for the 12 CO 7=3-2/ J=l-0 map of NGC 628, 
the H 2 (CO 7=3-2) gas mass fraction map for this 
galaxy is limited by the BIMA SONG 12 CO 7=1-0 
map. The THINGS moment maps for NGC 628 
show the Hi intensity dips in the central ~1' 
(~2.1 kpc) and there is little correlation between 
the 12 CO 7=3-2 and Hi in this area; beyond this 
radius they both trace the spiral structure (though 
this goes beyond the useful area of the 12 CO 7=1- 
data). In the limited regions where we have all 
three gas tracers we find the dense molecular gas 
is only the dominant gas phase at a few peaks in 
the map. It is difficult to draw any conclusions 
from this map given the limited overlap between 
the three gas tracers. 

In NGC 3521, the THINGS Hi map shows 
a larger central depression, that encompasses 
most of the region with molecular gas detections. 
Molecular gas in general dominates the inner re- 
gions of the galaxy (out to ~1.'5 from the center 
along the major axis, or about 4.7 kpc at the 
adopted distance). Dense/ warm molecular gas 
appears to be the dominant gas phase in several 
large patches of the gas disk, with a strong peak 
just south of the galaxy center where the dense 



14 



molecular gas is more than three quarters of the 
total gas along the line-of-sight, and a similar 
region further out to the south. The most no- 
table location with little dense molecular gas is 
a small region slightly north of the center, in ap- 
proximately the same location we see a dip in the 
12 CO J=3-2 distribution map (Fig. Eh). 

Molecular gas is in general the dominant form 
of gas in NGC 3627, in particular throughout the 
bar and most of the spiral arms. The 12 CO J=3- 
2 detections are particularly strong at four peaks; 
in the central starburst region, at both ends of 
the bar, and a peak to the south between one of 
the bar ends and the south-western arm. Atomic 
gas on the other hand is mostly absent in the bar 
except at the north-western end, but traces most 
of the spiral arms. As a result the map shows 
dense/warm gas is the dominant gas phase at the 
four peak locations (in the case of the central star- 
burst region overwhelmingly so), but not through- 
out the rest of the bar and spiral arms. The map 
strongly resembles the 12 CO J— 3-2/ J— 1-0 map 
for this galaxy, indicating that cooler and/or lower 
density molecular gas is the dominant phase out- 
side the four peaks. 

Fig. [7] shows the pixel-to-pixel variation in the 
12 CO J=3-2/J=l-0 ratio versus H 2 gas surface 
density (from 12 COJ=3-2, inclination and he- 
lium corrected) for all three galaxies. Through- 
out the vast majority of all three galaxies the ra- 
tio is less than 0.5. There is a discernible trend 
for NGC 3627, with higher surface density regions 
having higher ratios. This trend is less obvi- 
ous for NGC 3521, and not clear at all from the 
limited NGC 628 data. In general, while regions 
with low H2 (CO J=3-2) gas surface density can 
have a wide range in 12 CO J=3-2/J=l-0, regions 
of high H 2 (CO ,7=3-2) gas surface density only 
show high ratios. This fits with our assumption 
that 12 CO J=3-2 tends to trace warmer, denser 
regions than 12 CO J = l-0. Given that starburst re- 
gions in galaxies are often o bserved to have hig her 



12 CO J=3-2/J=l-0 ratios (|Meier et all 120011) . it 



is somewhat surprising that the starburst galaxy 
NGC 3627 generally has the lowest ratio of these 
three galaxies at a given molecular gas surface den- 
sity. 

Fig. [8] shows the pixel-to-pixel variation in 
the H 2 (CO J=3-2) gas surface density fraction, 
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Fig. 7. — Plot of the pixel-to-pixel variation 
in the 12 CO J=3-2/ J=l-0 ratio versus the dense 
H2 (CO J=3-2) gas surface density for all three 
galaxies. 12 CO J=3-2/ J=l-0 is calculated by di- 
viding our integrated 12 CO J=3-2 intensity map 
(main beam temperature scale) by the corre- 
sponding integrat ed 12 CO J=l-0 intensity (from 
Kuno et al.l (|2007l) . except for NGC 628 (*) where 
the BIMA SONG 12 COJ=1-0 data were used 
instead). S Hs ,( CO j=3-2) is calculated from our 
NGLS 12 CO J=3-2 observations using the as- 
sumption for dense gas ( 12 CO J=3-2/ J=l-0 = 0.6, 
see text), and has been corrected for helium frac- 
tion and inclination (using the inclina tions from 
Walter et alJl2008t Ide Blok et al.ll2008h . NGC 628 
is marked with solid black triangles, NGC 3521 
with open red circles, and NGC 3627 with open 
blue squares (these symbols are used for all subse- 
quent plots). We have only included pixels where 
the signal-to-noise for the NGLS 12 CO J=3-2 data 
is greater than 3. The error bars indicate a typical 
calibration uncertainty. 
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Fig. 8. — Plot of the pixel-to-pixel variation 
in the H2 (CO J=3-2) gas surface density frac- 
tion (S H2 (co J=3-2) divided by S H i+h 2 (co j=i-o)) 
verses total gas surface density (Hi mass + H2 
[CO J= 1-0] mass ) for t he two galaxies where we 
have iKuno et al.l (|2007t l 12 CO J =1-0 maps avail- 
able. Symbols are as in Fig.[7J All gas surface den- 
sities were corrected for helium content, and we 
applied an inclination correcti on where appropri- 
ate (using the inclin ations from lWalter et aO2008t 
de Blok et al.ll20081 ). We have only included pixels 
where the signal-to- noise for the NGLS 12 CO J=3- 
2 data is greater than 3. The error bars indicate 
a typical calibration uncertainty. 



tal gas surface density £hi+h 2 (CO j=i-o) (inclina- 
tion and helium corrected) f or the two galaxies 
where we have IKuno et al.l (|2007l ) 12 COJ=1-0 
maps available. Only pixels where the signal- 
to-noise for the NGLS 12 CO J=3-2 data is greater 
than 3 were included. The absence of points below 
~10 .M©pc~ 2 is most likely a result of the detec- 
tion limits of the three surveys, the signal-to-noise 
cutoff that we subsequently applied, and being 
biased towards pixels where 12 CO J=3-2 was de- 
tected in the NGLS. With the large scatter in the 
points it is difficult to discern any trend between 
the global gas surface density and the fraction of 
dense/ warm molecular gas, and in general the plot 
suggests there is no relation between gas surface 
density and gas compositi on. This may be c onsis- 
tent with results such as iLerov et al. (2008) with 
regards to the influence general ISM gas densities 
have on denser structures more closely associated 
with star formation, provided we can confirm it 
with measurements from other galaxies at a later 
stage. 

5. Gas Dynamics 

We compared the NGLS 12 CO J=3-2 velocity 
fields and velocity dispersion maps for our three 
galaxies to the THINGS Hi velocity fields a nd ve- 
locity dispersion maps (jWalter et al.l 120081) . The 
THINGS H 1 data cubes were first convolved with 
a Gaussian beam to the resolution of the JCMT 
12 CO J=3-2 maps (14'.'5), and then collapsed into 
moment maps, as described in § [2j To compare 
velocity fields we simply subtracted the NGLS 
12 CO J=3-2 fields from the THINGS Hi fields. To 
compare velocity dispersions we divided the NGLS 
12 CO J=3-2 velocity dispersion maps by the Hi 
velocity dispersion maps to get the ratio of the 
velocity dispersions. The resulting velocity field 
difference and velocity dispersion ratio maps are 
shown in the lower two panels of Figs. |U and[B] 
It should be noted that our linear resolution ele- 
ments for the three galaxies are of the order of sev- 
eral hundred parsecs, relatively large compared to 
the internal structures of the galaxies, and there- 
fore the gas traced by Hi and 12 CO J=3-2 at the 
same location in the maps may be physically some 
distance apart. 

Over the regions in NGC 628 where we detect 
12 CO J=3-2, the velocity field differences are gen- 
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erally small (absolute differences almost entirely 
< 6 kms -1 ) in comparison to the other two galax- 
ies, no doubt in part due to the low inclination of 
this galaxy (and hence small component of the ro- 
tation along the line of sight). There does appear 
to be a slight gradient across the field, with the 
northern (receding) end of the galaxy having H I 
velocities marginally greater than 12 CO J=3-2 ve- 
locities, and the southern (approaching) end show- 
ing slightly lower Hi velocities than 12 CO J— 3-2 
velocities. The extent to which this gradient is 
real is not clear, and would require a more exten- 
sive tracer of the molecular gas to expand the area 
where we can compare the velocities over the very 
extended Hi disk. The 12 CO J=3-2 velocity dis- 
persion is significantly below that of the atomic 
gas throughout the entire region where we detect 
12 CO J=3-2, with few regions above 50% of the H I 
velocity dispersion (the median ratio of the disper- 
sions is 34%). In general we expect the molecular 
gas to have lower velocity dispersions because it is 
confined to a colder, thinner disk. 

There are large regions of disagreement between 
the atomic and molecular velocity fields for the 
highly inclined galaxy NGC3521, with peak dif- 
ferences of up to ±40 kms -1 in the regions to the 
north and south of the galaxy center. Unlike the 
gradient seen for NGC 628, the differences appear 
to be more clump-like, and are in the opposite 
sense to that galaxy (higher 12 CO J=3-2 veloci- 
ties on the receding side, and vice versa). Exami- 
nation of a position-velocity slice along the major 
axis through the THINGS H I cube for NGC 3521 
show significant quantities of low velocity atomic 
gas. This gas following a non-circular rotation 
pattern, which is not detected in the 12 CO J=3-2 
data cube, appears to account for the large dif- 
ferences in the Hi an d 12 CO J=3-2 velocity fields. 
de Blok et al. ( 20081 and references therein) dis- 
cuss removal of such low velocity gas from the 
moment maps for THINGS, so that velocity fields 
best represent the rotational pattern of the galaxy. 
Presumably, if we compared an H I velocity field 
that better reflected the circular rotation with our 
12 CO J=3-2 velocity field, there would be fewer 
large differences. Fig. [5ji shows the 12 CO J=3- 
2 to Hi velocity dispersion ratios for NGC 3521. 
The 12 CO J=3-2 velocity dispersion is again much 
lower than for the atomic gas, of similar order to 
NGC 628. The median ratio of the velocity disper- 



sions is 30%, though again these velocity disper- 
sions ratios are likely affected by the non-circular 
H I gas and may be more similar if only H I in cir- 
cular rotation were used. 

Fig. [6J: shows the velocity field differences for 
NGC 3627. The molecular and atomic velocity 
fields for this galaxy are well matched, with the 
exception of the region around the nucleus and the 
southern end of the bar where the 12 CO J=3-2 ve- 
locities are over 30 kms -1 higher in some parts. 
This is the receding side of the galaxy, though 
it is notable that this region has relatively little 
Hi emission. The bar would also strongly affect 
the dynamics in this region. Fig. [6ji shows the 
12 CO J=3-2 to Hi velocity dispersion ratios for 
NGC 3627. The 12 CO J=3-2 velocity dispersion is 
closer to the H I velocity dispersions for this galaxy 
compared to the others, though most of the galaxy 
still has 12 CO J=3-2 velocity dispersions < 75% 
of the H I values. The median ratio of the veloc- 
ity dispersions is 42%. Most notable is the nu- 
cleus where the 12 CO J=3-2 velocity dispersion is 
in fact greater than Hi values. As noted in § |4] 
the molecular gas dominates the central starburst 
region, making up almost all of the gas (95% to 
99%) at the nucleus, so it is likely that the atomic 
velocity dispersion here is based on low signal-to- 
noise data. 

6. Star Formation 

6.1. Star Formation Rate and Gas Deple- 
tion Time 

To make a comparison between the present star 
formation rate (SFR) and our measurements of 
dense/warm molecular gas derived from NGLS 
12 CO J=3-2 observations we have produced maps 
of the rate of star formation from the available an- 
cillary data. We followed the same procedures as 
Paper I to produce SFR maps, and maps of the gas 
depletion time for H2 (CO J=3-2). Star formation 
rates were derived from the available SINGS Ha 
and 24 ixm m a ps us ing the procedure derived by 



Calzetti et all (|2007l . eqn. 7) 



SFR^oyr" 1 ) = 5.3 x 10" 42 [L(Ha) obs + 
(0.031 ±0.G06)£(24/zm)], 



(1) 



where the Ha and 24 /im. are both in erg s 1 , and 
L(24^m) is expressed as vL{y). This calculation 
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assumes an initial mass function (IMF) consisting 
of two power laws, unlik e the similar derivation in 
iKennicutt et al. (1998a) that uses a Salpeter IMF. 
By using both Ha and 24 /zm data we account 
for both unobscured and obscured star formation 
within the galaxies. 

The gas depletion time, which is also the inverse 
of the star formation efficiency (SFE), is simply 
the mass of the gas in question divided by star for- 
mation rate. This assumes that the gas involved 
will eventually be used up in the process of forming 
stars, and that no other process is replenishing this 
gas. While these may hold somewhat for the total 
gas population over the lifetime of the galaxy (if 
there is no external influence), they are unlikely to 
be the case for the dense/warm molecular gas we 
trace with 12 CO J=3-2. If we are only consider- 
ing how quickly the current supply of gas is being 
consumed for star formation, and not the longer 
term implications of when that gas might run out, 
this should not be a concern. Using our calcula- 
tion for the dense/ warm molecular gas mass from 
the NGLS 12 CO J=3-2 maps (H 2 [CO J=3-2]), we 
have produce maps of the depletion time of the gas 
presumably closest to the star formation. 

The maps of SFR and H 2 (CO J=3-2) gas de- 
pletion time for NGC 628 are shown in Fig. [9] Star 
formation in this galaxy appears to be mostly con- 
centrated around the spiral arm structure. In the 
center of the galaxy the SFR is generally lower 
compared to the sites of strong star formation in 
the surrounding arms. Accordingly this region 
has a relatively long gas depletion time (>1 Gyr) 
in comparison to much of the rest of the galaxy, 
though some sites along the arms have similar de- 
pletion times. 

Fig. M shows the SFR and H 2 (CO J=3-2) gas 
depletion time maps for NGC 3521. Star forma- 
tion appears to be strongly concentrated in the 
center of this galaxy, though there are likely strong 
projection effects in this near edge-on galaxy. The 
strongest SFR peaks are in the inner ~1' or so 
(~3 kpc), the same region in the inner parts 
of the 12 CO J=3-2 intensity map that shows a 
slight depression. This distribution is reflected 
in the gas depletion time map, where the deple- 
tion time in the center is less than 0.5 Gyr in 
some places, whereas some surrounding regions 
have longer times of 1-1.5 Gyr. Otherwise the de- 
pletion times for dense H2 gas in NGC 3521 are 
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Fig. 9. — a) Star formation rate map for NGC 628, 
derived from the SING S Ha and 24 /mi images 
(jKennicutt et alj|2003h . convolved to the JCMT 



beam at 850 /mi. Contour levels are 0.002, 0.004, 
0.008, 0.016, 0.032, and 0.064 Mgyr^kpc -2 . b) 
H 2 (CO J=3-2) gas depletion time (gas mass di- 
vided by SFR) for NGC 628 in Gyr. The magenta, 
black, and cyan contours mark where the gas de- 
pletion time is 0.5, 1.0, and 2.0 Gyr, respectively. 
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Fig. 10. — a) Star formation rate map for 
NGC3521, de rived from the SING S Ha and 
24 /Ltm images (jKennicutt et al.ll2003l) . convolved 
to the JCMT beam at 850 /im. Contour lev- 
els are 0.002, 0.004, 0.008, 0.016, 0.032, and 
0.064 A^oyr-ikpc" 2 . b) H 2 (CO J=3-2) gas 
depletion time (gas mass divided by SFR) for 
NGC3521 in Gyr. The magenta, black, and cyan 
contours mark where the gas depletion time is 0.5, 
1.0, and 2.0 Gyr, respectively. 



relatively uniform at ^1.0 ±0.5 Gyr. 

Fig. CCD shows the SFR and H 2 (CO J=3-2) gas 
depletion time maps for NGC 3627. There is in- 
tense star formation activity at both ends of the 
bar in this galaxy, as well as in the center of 
the bar. Other sites of strong star formation are 
scattered along the spiral arms. What appears 
to be a star formation peak showing up on the 
north-western corner of the field, outside the re- 
gion where we see 12 CO J=3-2 and other gas trac- 
ers, is not related to NGC 3627 but is an artifact 
from an over-exposed foreground star on the Ha 
image. The gas depletion time map shows that 
throughout the spiral arms the depletion time is 
similar to the other two galaxies (<2 Gyr). How- 
ever, in the middle of the bar around the central 
starburst region, despite the strong star formation 
here, the depletion time rapidly lengthens due to 
the high surface density of H 2 (CO J=3-2) here, 
and is over 4 Gyrs in the center pixel. 

With the exception of regions like the center of 
NGC 3627, the typical H 2 (CO J=3-2) gas deple- 
tion times throughout all three galaxies is of the 
order of 1 Gyr. Median gas depletion times are 

1.1 Gyr for NGC 628, 1.3 Gyr for NGC 3521, and 

1.2 Gyr for NGC 3627 (mean values are 1.2 Gyr, 
1.8 Gyr, and 1.4 Gyr, respectively). These are sim- 
ilar to the mean global gas depletion times found 
for three of the Virgo spiral galaxies studied in Pa- 
per I (1.1 to 1.7 Gyr). As we discussed in that pa- 
per, these times are comparabl e to those fo u nd for 
the spiral galaxies stu died in Lerov et al. ( 2008f) 
and Bigiel et al] (J2008) , but much shorter than a 
Hubble time, though as already mentioned they 
should not be interpreted as the time that the fuel 
for star formation will be exhausted. The relative 
consistency of the global gas depletion times for 
dense/warm molecular gas between spiral galax- 
ies with very different structures, gas composition, 
and environment is likely more an indication of the 
transitional nature of such gas in spirals rather 
than the properties of the galaxy itself. 

6.2. Star Formation Efficiency Relation- 
ships 

Using the Hi data from THINGS and other 
ancillary data from related surveys (SINGS, 
12 CO J=2-l f rom HERACLES, and GALEX UV 
observations) . [Lerov et al" ( 2008f ) have performed 
an extensive analysis of the star formation effi- 
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Fig. 11. — a) Star formation rate map for 
NGC 36 27. derived from the SIN GS Ha and 24 ^m 

convolved to the 



2003) 



images (jKennicutt et al. 
JCMT beam at 850 ^m. Contour levels are 
0.002, 0.004, 0.008, 0.016, 0.032, 0.064, 0.128, 
and 0.256 .A/( yr" 1 kpc" 2 . b) H 2 (CO J=3-2) 
gas depletion time (gas mass divided by SFR) for 
NGC 3627 in Gyr. The magenta, black, cyan, and 
green contours mark where the gas depletion time 
is 0.5, 1.0, 2.0, and 4.0 Gyr, respectively. 



ciency within a sample of nearby galaxies. This 
included examinations of various star formation 
laws, such as the we l l know Schmidt-Kennicutt 
" 1998allbh . In Fig. [lj we 
two plots based on th e SFE ver- 
~" Fig. 5. 



law (jKennicutt et al 
have shown 



sus E gQS plots shown in iLerov et al.l Fig. 5. We 
have reproduced these plots for our galaxies on 
a pixcl-to-pix el basis (azimuth ally smoothed data 
are used in ILerov et al.1 120081 ). For reference we 



have included the same lines shown in ILerov et al 



in both panels, marking the Hi-to-H 2 transition 
surface density in spirals (vertical dotted line at 
^gas = 14 A^opc -2 , only relevant where Hi is in- 
cluded) , and the relation SFE oc E° a 5 s (an approx- 
imation to the Schmidt-Kennicutt star formation 
law). 

The left panel (Fig. [Hi) shows the H 2 (CO J=3- 
2) gas star formation efficiency (star formation 
rate per unit of dense H 2 gas) versus the H 2 
(CO J=3-2) gas surface density for all three galax- 
ies. The H 2 (CO J=3-2) gas surface density used 
in both values was estimated using the assumption 
for dense/warm molecular gas mentioned earlier 
( 12 CO J=3-2/ J=l-0 = 0.6), and hence we are only 
looking at the efficiency of star formation in the 
molecular gas regions most closely associated with 
star formation. The right panel (Fig. \T2b) on the 
other hand show the behavior of the more global 
gas population. It shows the total gas (Hi + H 2 
[CO J=l-0]) star formation efficiency versus to- 
tal gas surface density (£hi+h,(co j=i— o))- Hence 
much of the gas traced in this plot is much less 
directly linked to star formation than in the first 
panel (and does not use our NGLS 12 CO J=3-2 
data) , though th e plot is more direc tly comparable 
to the results of Lerov et alj (|2008T ) . 

Both plots have large scatter and it is difficult to 
discern clear trends in either. The left panel shows 
that the star formation efficiency of H 2 (CO J=3- 
2) gas is relatively high (as indicated in the previ- 
ous section by the short gas depletion times). The 
efficiency is relatively constant with gas density, 
and if anything becomes less efficient when the H 2 
(CO J=3-2) gas surface density is higher, and it 
does not follow the line of the Schmidt-Kennicutt 
law. However, given the order of magnitude peak- 
to-peak scatter of the distribution it is difficult to 
conclude whether the efficiency is constant or de- 
clining with gas density. In the right panel we see 
the star formation efficiency is about an order of 
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Fig. 12. — Pixel-to-pixel plots of star formation efficiency (inverse of the gas depletion time) versus gas 
surface density (symbols as in Fig. [7]) . a) The H2 (CO J— 3-2) gas star formation efficiency (star formation 
rate per unit of H2 [CO J=3-2] gas) versus H2 (CO J=3-2) gas surface density for all three galaxies, b) The 
total gas (Hi + H2 [CO J— 1-0]) star formation efficiency versus total gas surface density (£hi+h 2 (co j=i-o)) 
for the two galaxies where we have iKuno et al.1 fl2007T) 12 C O J =1-0 maps available. In both panels, the 
lines correspond to the same markings in Lerov et al.l (|2008 ) Fig. 5. The vertical dotted line marks S oas = 

2008). 



14±6A^opc iLerov et al 



14 A^opc , the Hi-to-H2 transition surface density in spirals (S 9 , 
The dashed line shows the relationship SFE oc S° a 5 s , an approximation to the Schmidt-Kennicutt star 
formation law. All gas surface densities were corrected for helium content, and we applied an inclination 



correction where appropriate (using the inclinations from Walter et al. 20081 de Blok et al. 2008f) . For both 



plots we have only included pixels where the signal-to- noise for the NGLS 12 CO J =3-2 data is greater than 
3. The error bars indicate a typical calibration uncertainty. 
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magnitude lower than the first plot, as might be 
expected when the total gas content is taken into 
account. The absence of points at lower densi- 
ties is likely due to the bias towards regions where 
12 CO 7=3-2 was detected. There is no clear trend 
in this plot due to the large scatter and narrow 
range of gas surface density. Overall there is little 
evidence that either plot is following a Schmidt- 
Kennicutt star formation law. 

Fig. ITBI shows a plot of SFE H2 ( CO j=3_ 2 ) versus 
the 12 CO 7=3-2/ 7=1-0 ratio for the two galax- 
ies where we have iKuno et al.l (|2007j ) 12 C0 7=1- 
maps available. The scatter is again large but 
maybe a bit tighter than the previous plots, and 
for the most part the SFE appears to be flat 
or slightly declining with increa sing 12 CO 7=3- 
2/J=l-0. iMuraoka et~aH (|2007t ) also looked at 
the relationship between SFEh 2 and 12 CO 7=3- 
2/J=l-0 for M83 (NGC5236), and found a clear 
trend for higher SFE at higher 12 CO 7=3-2/ 7=1- 
ratios. Their azimuthally smoothed data only 
cover a small range in 12 CO 7=3-2/7=1-0 com- 
pared to our pixel-to-pixel measurements, and 
their SFEh 2 range is close to our scatter in that 
value, but our results do not seem to support this 
trend. The most likely explanation from our plots 
is that the star formation efficiency in theses two 
galaxies is independent of the 12 CO 7=3-2/7=1-0 
ratio. 

Fig. ITU shows a plot of SFE H2 ( CO j=3-2) ver- 
sus the H 2 (CO 7=1-0) gas mass surface density 
fraction, E H , (r,n 3=3-2) /^HI+ H^CQ J=] -n) , for the 
two galaxies where we have IKuno et al.l (2007) 
12 CO 7=1-0 maps available. As with the previous 
plots there is a large scatter in the points close 
to a magnitude in star formation efficiency. But 
as with the previous plot, and the first panel of 
fig. |T2l as the surface density of the dense/ warm 
molecular gas increases the star formation effi- 
ciency remains constant or becomes less efficient. 
This could indicate that in regions of dense and / or 
warm molecular gas the efficiency is either inde- 
pendent of the gas density, or become slightly 
less efficient as the density increases. Investiga- 
tions with larger samples with extensive multi- 
wavelength data are needed to confirm if either 
is the case. 
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Fig. 13.— Plot of the H 2 (CO 7=3-2) gas star 
formation efficiency (star formation rate per unit 
of H 2 [CO 7=3-2] gas) versus 12 CO 7=3-2/7=1- 
for the two galaxies where we have IKuno et al 



(2003) 12 CO 7=1-0 maps available (symbols as in 
Fig. E]). 12 CO 7=3-2/7=1-0 is calculated by di- 
viding our integrated 12 CO 7=3-2 intensity map 
(main beam temperature scale) by the corre- 
sponding integra ted 12 CO 7=1-0 intensity (from 
IKuno et alJ[2007h . The H 2 (CO 7=3-2) gas sur- 
face density for the star formation efficiency cal- 
culation was estimated using the assumption for 
dense gas ( 12 CO 7=3-2/7=1-0 = 0.6, see text). 
All gas surface densities were corrected for he- 
lium content, and we applied an inclination cor- 
recti on where appropriate (using the inclina tions 



from lWalter et alj|2008t Ide Blok et ai] |2008). We 



have only included pixels where the signal-to-noise 
for the NGLS 12 CO 7=3-2 data is greater than 3. 
The error bars indicate a typical calibration un- 
certainty. 
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Fig. 14.— Plot of the H 2 (CO J=3-2) gas star 
formation efficiency (star formation rate per unit 
of H 2 [CO J=3-2] gas) versus the H 2 (CO J=3- 
2) gas surface density fraction (Eh 2 (co j=3-2) di- 
vided by Em+Ho fno ,i=i -ni) f or th e two galax- 
ies where we have lKuno et al.l (|2007l) 12 CO J=l-0 
maps available (symbols as in Fig. [7]). The H 2 
(CO J=3-2) gas surface density for the star for- 
mation efficiency calculation was estimated using 
the assumption for dense gas ( 12 CO J=3-2/J=l-0 
= 0.6). All gas surface densities were corrected for 
helium content, and we applied an inclination cor- 
recti on where appropriate (using the inclina tions 
from lWalter et all 120081; Ide Blok et al.ll2008l ). We 
have only included pixels where the signal-to-noise 
for the NGLS 12 CO J=3-2 data is greater than 3. 
The error bars indicate a typical calibration un- 
certainty. 



We have used the JCMT to map the 12 CO J=3- 
2 emission from three non-cluster spiral galaxies, 
completed as part of the ongoing Nearby Galax- 
ies Legacy Survey. T hese galaxies are all inc luded 
in the SINGS survey (|Kennicutt et al.ll2003l) . and 
as such have excellent published data sets cov- 
ering their dust and ISM properties. Combined 
with the new NGLS 12 CO J=3-2 data, these data 
set also allow us to probe a range of galaxy 
properties including differences in the distribu- 
tion of 12 CO J=3-2, 12 CO J=l-0, and Hi emis- 
sion, galaxy dynamics, and relationships between 
star formation rates and warm/dense molecular 
gas. 

These galaxies all have moderate to strong 
12 CO J=3-2 detections over large areas of the 
fields observed by the survey, showing resolved 
structure and dynamics in their warm/dense 
molecular gas disks. We find some large differ- 
ences between the 12 CO J=3-2 and Hi velocity 
fields of some galaxies, in particular NGC3521. 
These differences appear to be due to the pres- 
ence of large quantities of low velocity H I gas in 
the THINGS data cubes. 

Global gas depletion times for dense/ warm 
molecular gas in these galaxies, as well as those 
examined in Paper I and other surveys, are rela- 
tively consistent with each other (in the range 1 to 
2 Gyr). This is despite very different structures, 
gas compositions, and environments. This is likely 
more an indication of the transitional nature of 
such gas in spirals rather than the properties of 
the galaxy itself. 



Similar to the results from the iLerov et al 



(|2008l ), we do not see any correlation of the star 
formation efficiency with the gas surface density, 
with at best a slight decline in efficiency in high 
gas density regions traced by 12 CO J=3-2. Ad- 
ditionally, we find that the star formation effi- 
ciency of the dense/warm molecular gas traced 
by 12 CO J=3-2 is flat or slightly declining as 
the 12 CO J=3-2/ J=l-0 ratio increases, in con- 
trast to the c o rrelat ion found in the study by 
Muraoka et ail (|2007l ) into the starburst galaxy 
M83. A similar weak trend is also seen when ex- 
amining star formation efficiency versus the H 2 
(CO J=l-0) gas mass surface density fraction. All 
three of these could indicate that in regions of 
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dense and/or warm molecular gas the efficiency is 
either independent of the gas density, or become 
slightly less efficient as the density increases. 
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Fig. 15. — The basic layout of the sixteen recep- 
tors in the HARP-B array. Each receptor sits in 
the center of a square. The layout is shown tilted 
by 14?04 clockwise, one of the angles relative to 
the scan axis used for scanning. Depending on the 
part of the sky in which observations are taken, it 
may also be orientated 14?04 the other way, and 
the use of a K mirror in the telescope introduces 
further orientation flips of 90°. The two grey re- 
ceptors (H03 and H14) were missing or unusable 
during all observation. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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A. NGLS Raster Map Observations and their Reduction Strategy 

This section describes the observing and data reduction methods used for the NGLS HARP-B 12 CO J=3-2 
raster observations of NGC 628, NGC 3521, and NGC 3627. This description is also relevant for the reduction 
of other raster map data taken for the NGLS, such as the Virgo spiral galaxies in Paper I (NGC 4254, 
NGC 4321, NGC 4569, and NGC 4579). It can also be generalized for jiggle map data observations (which 
often require less or no flagging). 

After a series of Science Verification runs between April and November 2007, the N GLS started takin g 



science data using the 16 element HARP-B array with the backend spectrometer ACSIS ( Smith et alj |2003). 
Since that time there has been approximately one dedicated JLS run per month that the three surveys share. 
Observations for the three galaxies we are examining here, NGC 628, NGC 3521, and NGC 3627, were taken 
over multiple JLS runs between November 2007 and March 2008. 

The James Clerk Maxwell Telescope (JCMT) is a 15 meter single dish submillimeter telescope on Mauna 
Kea, Hawai'i. It is currently the largest single dish submillimeter telescope in the world. HARP-B is a 
16-element sparsely filled array of heterodyne receptors operating in the submillimeter window centered at 
850 /im (350 GHz). The receptors are arranged in a 4x4 pattern and separated by 30", approximately two 
beam- widths (~14'/5 at this wavelength). See Fig [15] for the basic layout of the array (with two receptors 
that were missing during observing marked). To produce a fully sampled map we need to set up a Minimum 
Schedulable Block (MSB) that either scans or jiggles the array around the field. For all three galaxies in 
this paper, which have fields larger than the array footprint of 2', we used a raster scanning technique for 
mapping. 

The default HARP-B raster scanning mode is to rotate the array by ±14?04 to the scan axis and to scan 
along the long axis of the field, which produces scan rows at 7'.'2761 intervals. On the outer edges of the 
field coverage is sparse where half the array receptors start already within the target field, so we added an 
additional 60" to the start and end along the direction of the scan to account for this. There are significant 
variations in signal-to-noise between receptors in the array, which can produce striping in raster maps where 
different rows have different sensitivity. To reduce this effect and produce maps with more even sensitivity 
coverage we used a 'basket weave' technique, scanning along both the major and minor axes of the field. 
Additionally, we found in the science verification phase of the survey that missing and bad receptors in the 
array caused large variations in our completeness across the target field, and even holes in our coverage 
where missing/flagged scans crossed in the basket weave. To counter this problem, we used a half-step scan 
method, only moving half an array step (58'.'2) between scan rows/columns. In this way, each position was 
covered by two receptors in both scan directions. However, using the half-step method does require that an 
additional 60" be added at either end perpendicular to the scan direction so that all points within the target 
field have the same integration time. With the previous addition of 60" along the scan axis, we are adding 
60" around all four sides of the field for both scan directions. 

Each galaxy's MSB was set up with a sample time per position of 2.5 seconds, so with the half-step basket 
weave pattern each pointing will have 10 seconds integration time per MSB (if none of the receptors that 
contribute to that pointing are flagged or missing). We repeated these MSBs up to seven times until we 
reached the survey's target RMS (19 mK after binning to 20 kms -1 resolution). The backend spectrometer 
ACSIS was set up with a bandwidth of ~1 GHz and a resolution of 488.28 kHz (0.423 kms -1 at the frequency 
of the observations), centered at the systemic velocity for the galaxy being observed (where the 12 CO J=3- 
2 line rest frequency is 345.79 GHz). Appropriate pointing and calibration observations were taken during 
observing nights. Individual setup and observing details for the three galaxies in this paper are in summarized 
Table El 

Data reduction and most of the analysis were carried out with recent JAC maintained versions of the 
Starlink software collection (initially the Humu version, and later Lehuakona). Most of the Starlink tasks 
we used in the reduction were in the KAPPA (most reduction and analysis tasks) and SMURF (cube 
making) packages, with some tasks in CUPID (clump finding) and CONVERT (file exporting) also used. 
The visualization programs GAIA and SPLAT were used to inspect the raw and processed data. Additionally, 
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analysis tasks in MIRIAD were used after the creation of cubes and moment maps. The general reduction 
procedure was to inspect the raw data files, flag poor data where necessary, combine the raw data into a 
cube, subtract a baseline from the cube, and collapse it into moment maps. Details of these steps follow. 

Throughout the period that the survey observations were taken, the state of the 16 receptors in the array 
varied. Two receptors (H03 and H14) were unavailable for the duration of the project and for the most part 
were automatically flagged at the telescope, while several other receptors had various, mostly intermittent 
problems that led to bad baselines, spikes, and false lines. We used the task chpix (KAPPA) to flag such 
receptors where required. Additionally, almost all scans contained a narrow interference spike that was also 
flagged. This spike occurred at the same frequency on all receptors (with varying strength) throughout 
a scan, and would remain at that frequency over the course of a night or even over several runs. Some 
scans taken for NGC3521 were excluded when pointing observations taken straight after showed significant 
pointing errors (>5"), while all data for NGC628 taken on 11 January 2008 were excluded due to poor 
baselines and high noise on virtually all receptors. Other observations that did not meet the survey quality 
assurance criteria (such as high T sys or many receptors with bad baselines) were included, but often with 
heavy flagging of poor receptors. 

After flagging, the raw data were combined into cubes using the task makecube (SMURF). We produced 
cubes with 7'.'2761 pixels, and trimmed the velocity axis to the central ^800 kins" 1 of the bandpass. We 
used a sinc(7rx)sinc(fc7ra:) ('SincSinc') kernel in makecube for the pixel weighting function. The resulting 
cube was then trimmed to remove the leading and trailing scan ends outside the target field where our scan 
coverage is incomplete. After removing a simple first order polynomial baseline (fit to the line-free ends of 
the bandpass) and binning to ^20 kms -1 , this cube was examined for further bad baselines and to check if 
we had reached our target depth for the survey. If necessary the process of flagging and cube making was 
repeated until we were satisfied that the poor baselines had been removed. 

To produce a more accurate fit to the baseline level, we used a method in which we masked out emission 
in the cubes, and then fit third order polynomial baselines to the unmasked data. To make the baseline 
mask, we first smoothed the cube with a [3 pixel, 3 pixel, 25 channel] box filter, and then used the task 
mfittrend (KAPPA) to make a mask of all of the features in the smoothed cube. In automatic mode, the 
'method' parameter in mfittrend allows the choice of several masking methods that have varying success at 
picking up broad and narrow line features. As the galaxies we are looking at contain a mixture of narrow 
and broad 12 CO J=3-2 lines, it was important to produce a mask that contains all emission features. 

Individual methods can miss some emission, such as the default 'Single' method that picks up narrow lines 
well, but can often miss broad features. To produce a better mask, we combined the results of two of these 
methods, 'Single' and 'Global' (which does a better job on broad lines, but picks up more noise especially 
at the line- free ends of the bandpass). We altered this combined mask to leave the line- free regions at the 
ends of the bandpass unmasked (determined by considering the published width of the H I line and visual 
examination of the cube). Additionally, for NGC 3521 and NGC 3627, which have broad lines in their centers, 
we included a block mask over a small central region for the full width of the line so that only the line-free 
ends would be used for the baseline fit. The block mask for NGC 3521 was split into overlapping northern 
and southern components that cover different velocity ranges. The final combined masks were applied to the 
original cube, and a third-order polynomial was fit to the masked cube to produce the baseline fit. These 
baselines were subtracted from the original cube to make the final baseline subtracted cube. 

Moment maps of the galaxies were made from the base line subtracted cubes using the Starlink imple- 
mentation of the clumpfind algorithm ( Williams et al.lfl994h . the Cupid task findclump, to identify emission 
above 3a. Sensitivity varied across the field, particularly in the regions outside the target field (see above), 
so before running it through the moment map procedure we first scaled the cubes by the noise level in order 
to better identify real clumps. To do this, we produced a map of the standard deviation in the line-free 
channels of the cube (using the baseline mask we previously made), and divided the cube by this map. This 
scaling was removed again after producing the maps so that the units were correct. 

As with the baseline mask steps, we made a smoothed cube with a [3,3,25] box filter. Using the baseline 
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mask previously made for the baseline fit, we calculated the RMS in the line free regions of the cube. We 
then used findclump (Cupid) on the smoothed cube to find all emission greater than 3c, and create a mask to 
exclude the regions without emission. The clump mask was then applied to the original baseline subtracted 
cube, and the result was collapsed into two dimensional moment maps using collapse (KAPPA). We produced 
maps of the integrated intensity (moment 0), the intensity weighted velocity co-ordinate (moment 1, the 
velocity held), and the intensity weighted dispersion (moment 2, the velocity dispersion). We converted the 
12 CO J=3-2 data from corrected antenna temperature (T^) to the main beam temperature scale by dividing 
by f]MB = 0.6. This correction is different from the value of 0.67 used in Paper I, and all 12 CO J=3-2 
measurements from that paper used here are corrected for this difference. The resulting moment maps are 
discussed in § [3J 
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